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Abstract. The radio, optical, X-ray and gamma-ray nebulae that surround many 
pulsars are thought to arise from synchrotron and inverse Compton emission. The 
energy powering this emission, as well as the magnetic fields and relativistic particles, 
are supplied by a "wind" driven by the central object. The inner parts of the wind can 
be described using the equations of MHD, but these break down in the outer parts, 
when the density of charge carriers drops below a critical value. This paper reviews the 
wave properties of the inner part (striped wind), and uses a relativistic two- fluid model 
(cold electrons and positrons) to re-examine the nonlinear electromagnetic modes that 
propagate in the outer parts. It is shown that in a radial wind, two solutions exist 
for circularly polarised electromagnetic modes. At large distances one of them turns 
into a freely expanding flow containing a vacuum wave, whereas the other decelerates, 
corresponding to a confined flow. 
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1. Introduction 

The size and quality of the observational database relating to pulsar wind nebula? 
has increased dramatically in the past five years, stimulating research into both the 
modelling of their spectra and the theory of the winds themselves (for recent reviews 
see [HE])- The basic idea of the pulsar wind as a mix of particles, fields and waves, 
was outlined over thirty years ago [3]. Close to the pulsar, the wind moves outwards 
at relativistic speed. When the the ram-pressure drops to approximately that of 
the surroundings, a "termination shock" occurs, and strong synchrotron and inverse 
Compton radiation is emitted in the relatively slowly moving plasma in the surrounding 
"nebula". However, although progress has been made on understanding the structure 
of the nebula, jH El |6] there are still difficulties associated with the propagation of the 
wind. For example, the questions of where and how the Poynting flux that it carries is 
dissipated into particles remains open [7J El [101 [TT]. This paper, does not present a 
solution to these problems, but concerns itself with an analysis of two simple, nonlinear 
waves that propagate in the wind, and are likely to play a crucial role in shaping the 
global solution. 

A rotating magnetised neutron star can be expected to emit predominantly linearly 
polarised waves in the equatorial plane, and circularly polarised waves along the rotation 
axis. In addition, the waves in the equatorial region may carry a non-vanishing phase- 
averaged component of the magnetic field component. At intermediate latitudes, the 
waves will, in general, be elliptically polarised and may also carry a non-oscillatory 
component. A linear analysis of these waves is clearly inadequate, since the characteristic 
quiver frequency of an electron at points interior to the termination shock front is many 
orders of magnitude larger than the rotation frequency of the star. A nonlinear analysis, 
however, is feasible only for a few simple modes. 

The first mode we select for analysis is a linearly polarised wave — the striped wind 
- which has vanishing phase-averaged magnetic field. This wave can be thought of 
as an entropy wave or, alternatively, a rotational discontinuity. It requires a minimum 
plasma density in order to propagate, and so can exist only up to some maximum 
distance from the pulsar. Its subluminal phase velocity equals its group velocity, and 
the field is essentially magnetostatic in the wave frame. Since it propagates at highly 
super-magnetosonic velocity, its evolution is fixed by boundary conditions close to the 
star. This wave is likely to provide a good description of the interior zone of a pulsar 
wind in the equatorial regions. 

Following this, the propagation of nonlinear electromagnetic waves of superluminal 
phase velocity is addressed. These waves cut off if the plasma density exceeds a certain 
limit, and so propagate only beyond some minimum distance from the star. Both 
linearly and circularly polarised nonlinear plane waves can be constructed analytically, 
but the formalism is much simpler for the case of circular polarisation. In this paper, a 
new analysis is presented of circularly polarised waves in spherical geometry, using the 
short wavelength approximation. This is appropriate in the polar regions of a pulsar 
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wind. However, one can expect that the linearly polarised electromagnetic modes in 
the equatorial region exhibit qualitatively similar behaviour. Because the phase speeds 
are superluminal, both inner and outer boundary conditions are required to specify the 
wave structure. At large radius, two types of solution are found: in one of them the flow 
velocity stagnates in a region of constant density, whereas, in the other, it asymptotes to 
constant radial velocity. The former applies to pulsar wind nebulae, which are confined 
by an external medium. The latter applies to unconfined winds. 

2. Parameters of the wave 

Consider radially propagating waves in a region far outside the light-cylinder of the 
pulsar, i.e., f = r/r^ ^> 1, where tl = cP/2tt is the radius of the light cylinder and P is 
the pulsar period. In this region, the magnetic field is approximately transverse, because 
the longitudinal (radial) component decays as f ~ 2 , whereas transverse component decays 
as f _1 . Furthermore, a plane wave approximation is appropriate, since the wavelength is 
small compared to the curvature of the wavefront. Though approximately radial, pulsar 
winds are definitely not spherically symmetric, which is accounted for by specifying a 
finite solid angle Q s < Air in which the wind propagates. 

The wave properties depend on the phase-averaged radial particle flux density (J), 
the phase-averaged total radial energy flux density (F), (which is the component T 0,1 
of the sum of the stress-energy tensors of the fields and particles), and the phase- 
averaged flux density of radial momentum, (Q), (which is the component T 1 ' 1 ). This 
latter quantity depends on both (F) and on the way in which the energy flux is shared 
between the particles and fields. For example, for a vacuum wave, the momentum flux 
density equals the energy flux density/c, whereas for cold particles of speed v in the 
absence of fields the momentum density equals v/c 2 times the energy flux density. 

Instead of the fluxes, however, it is more convenient to use dimensionless 
parameters. Conventionally, a steady, relativistic MHD wind is characterised by the 
parameter \x introduced by Michel [12], which is the luminosity carried by the wind per 
unit rest mass. It gives the Lorentz factor each particle would have if all the energy flux 
were carried by the particles: 



In a pulsar wind, fi is constant, (provided pairs are not created or annihilated and no 
energy is dissipated into radiation in the region under consideration), but its value is 
uncertain. For the Crab Nebula the average of /x over the pulsar lifetime can be found 
by combining an estimate of the number of radiating electrons/positrons, obtained from 
the synchrotron luminosity, with an estimate of the average magnetic field, obtained 
from the gamma-ray luminosity. But the current value is thought by many to be rather 
higher than in the past. One estimate is /i ~ 10 4 [2] but the real value might be larger. 

The distribution of energy flux between particles and fields can be described by 
the magnetisation parameter a, defined, in an MHD flow, as the ratio of the enthalpy 
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density of the magnetic field to that of the plasma, measured in the plasma rest frame. 
This quantity depends on the fraction of the pulsar luminosity that goes into creating 
pairs, i.e., not only on their number, but also on the Lorentz factor with which they 
are injected into the wind, a is the ratio of the maximum possible Lorentz factor (the 
parameter /i) to this injection Lorentz factor. This is generally thought to be high 
~ 10 2 -10 3 , so that a w l(T 2 //-l(ry 

Finally, as a dimensionless expression of the total energy flux it is convenient to 
choose the strength parameter of a circularly polarised vacuum wave that carries the 
same energy flux as the nonlinear wave. Writing E e g for the magnitude of the electric 
field in this wave, the strength parameter is defined as a = eE e s/(maj), and the required 
value of E eS is determined from cE 2 s /(4:ii) = (F): 



47re 2 (F) 
m 2 c 3 u 2 
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In a pulsar, this quantity varies with position in the wind: a 
strength parameter evaluated at r = r L : 

1/2 
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a L /r where a^ is the 
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(3) 

where L is the luminosity carried by the (radial) wind, and L 38 = L/(10 38 ergs -1 ). For 
the Crab, L = 4.6 x 10 38 ergs _1 ; many pulsars have L = 10 31 — 10 33 ergs _1 , and some 
go down to 10 28 ergs _1 . The value of Q s is not known, but is often assumed to be ~ 1. 
In a steady, radial wind, 

r 2 (J) =Cj 



f 2 (F) =mc 2 C F 
f 2 (Q) = mcC Q 



(4) 
(5) 
(6) 



where Cj^q are constants and /i = CpjCj. The particle density is frequently expressed 
in units of the "Goldreich- Julian" density at the light cylinder [13]. For a magnetically 
dominated, relativistic flow, the corresponding "multiplicity" parameter is k = aL/(4/x). 



3. The striped wind — an MHD wave 

Coroniti [H] and Michel [15] suggested idealised solutions in which the wind settles down 
into stripes of magnetic flux of alternating polarity containing cold plasma, separated 
by current sheets containing hot plasma. In this solution, the plasma is at rest in a 
reference frame that, with respect to the star, moves radially outwards with Lorentz 
factor 7mhd = (1 — Pmub) ■ m the comoving frame, the electric field vanishes, and 
the magnetic field is toroidal. There are two current sheets per wavelength. 

Such solutions can exist only if the wind carries a sufficient number of particles 
to carry the current that Ampere's law requires must flow in the sheets. To within a 
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factor of the order of unity, this condition is equivalent to the requirement that particles 
can be confined in the sheet, i.e., that the gyro-radius of a hot particle moving in the 
magnetic field present in the cold part of the flow is smaller than the sheet thickness. 
It was noticed in [13] that this condition implies a minimum dissipation rate: as the 
density decreases radially, the minimum sheet thickness increases. According to [13J, 
dissipation forces the wind to accelerate, which, because of time-dilation, reduces the 
apparent dissipation rate observed in the lab. frame. Of course, the actual dissipation 
rate is unknown. If it lies below the minimum, then the solution moves out of the range 
of validity of the MHD equations. If it lies above the minimum, the sheet expands, 
consuming the surrounding cold plasma and annihilating the field embedded in it. In 
the equatorial plane, adjacent sheets carry equal and opposite magnetic fluxes, and it 
is conceivable that the wave eventually annihilates the entire magnetic field. However, 
at higher latitudes, a residual flux remains after removal of the alternating component. 

Several possible dissipation or "reconnection" scenarios were investigated in [9]. 
The conclusion, in the case of the Crab, is that the maximum possible dissipation 
rate, combined with a very high particle content (small /i), could marginally enable the 
conversion of all the wind energy flux into particle kinetic energy before the termination 
shock is encountered. Given a reconnection prescription, self-similar solutions can 
be found that describe how the flow accelerates. For example, for dissipation at the 
minimum rate, the radial dependence can be approximated by [9] 

7mhd = fi\JSr/(7ra L ) (7) 

This relation holds in the supersonic part of the flow (7mhd > M 1//3 ), and before the 
magnetic field is annihilated completely (7mhd </•*)• In terms of the radius, these 
restrictions imply > f > a^/jj 4 ^ 3 . 

However, there exists an additional, independent constraint that is imposed by the 
use of a short-wavelength approximation in this analysis. This requires that the striped 
wind pass through a series of quasi-stationary states as it expands in the radial geometry. 
For this to remain valid, the gyro-period of particles contained in the sheets must be 
short compared to the expansion timescale. An equivalent formulation of this condition 
is that the speed at which the sheet expands must remain subsonic. This implies an 
upper limit on the wind Lorentz factor, and, hence, the radius: 

7mhd < 7m ~ a L //x r < a\/ '// (8) 

above which the wave can no longer be described by the MHD equations. 

In order to relate the striped wind solutions, which apply in the inner part of the 
wind, to the electromagnetic waves to be described in the following section, which apply 
in the outer part of the wind, it is useful to write down the phase-averaged fluxes carried 
by the wind, assuming the contribution of the current sheets to be negligible. Noting 
that the electric and magnetic fields are related by E = (3 M n D B, are orthogonal, and are 
constant in between the current sheets, (although they change sign across each sheet), 
and that the wave speed equals the fluid speed, i.e., 7mhd = 7 = J 1 +P|| 2 ) where 7 is 
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the particle Lorentz factor and p\\ the radial momentum in units of mc, one finds: 

(J) = 2np\\c (9) 

B 2 

(F) = 2nmc 3 p in + c(3 MUD — (10) 

(Q) = 2nmcV + (l + /^hd) (11) 

where n is the proper number density of electrons, which equals that of positrons. In 
terms of the magnetisation parameter, a = B 2 / (87r7^ HD nmc 2 ): 

{l) = A* = 7M HD (l + a) (12) 



(13) 



mc 2 (J) Cj 

and so in this case, in which the effect of the current sheets is neglected, both 7mhd and 
u remain constant in the flow. For /i>l and o <C fi 2 , the energy and momentum flux 
densities are almost equal: 

C F Cq 1 + cr 

irr K cj""-^r (14) 

4. Electromagnetic waves 

A relatively simple description of the outer parts of a pulsar wind can be based on a 
model in which two cold, relativistic, fluids (electrons and positrons) propagate radially 
outwards in transverse electromagnetic fields. Strong waves in systems similar to this 
were analysed by [HI (TTJ [HI [19], using the method of [20]. In connection with pulsars 
[2Tj and |22j investigated strong, linearly polarised waves in spherical geometry. Here the 
effect of spherical geometry on circularly polarised electromagnetic waves is analysed. 

Under the above assumptions, the proper number density is the same for each fluid: 
n + = ri- = n. In spherical polar coordinates (r, 6, <fi) their radial velocity is also the 
same: v r+ = t> r _ = v r and the toroidal and azimuthal velocities are equal in magnitude 
but of opposite sign: vg + = —v e _, = — It proves convenient to use the Lorentz 
factor of the fluids 7+ = 7_ = 7 = cj ^Jc 2 — v 2 — v 2 + — v 2 + , to define the dimensionless 
radial momentum: p» = Vrj/c, and to introduce complex quantities representing the 
dimensionless transverse momentum: p±_ = {vg + + iv^ + ) 7/c (= — (vg_ + iv^J) 7/c) and 
the electric and magnetic fields: E = Eg + iE^, B = Bg + iB^. Then the equations to 
be solved are the continuity equation 
d , c d ( 2 \ 

dt^ + r^d^y m ) =0 (15) 

the fluid equations of motion 

^ = -—lm(p ± B*) (16) 
dr mc 

^ = — (jE + i Pll B) (17) 
dr mc v ' 



<j) = OJ 



t - 



(21) 
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P- =—Re( P± E*) (18) 
dr mc 

(where r is the proper time (dt = 7<ir) and d/dr = d/dr + cp\\d/dr is the convective 
derivative), together with Faraday's law and Ampere's law: 

1 d . „. i dB , , 

( r #_ =0 19 
r ar c at 

If) i r)W 

-^-(rB) + -^- = -t8nnep ± (20) 
r or c at 

In the limit of large r, approximate solutions to these equations can be found that 
describe nonlinear waves, using the method of [20J. Introducing the phase 

1 r r dr' 
cJo P{r>)_ 

one can identify a circularly polarised solution valid for superluminal phase speeds 
((3 > 1). In it, the fields are 

B = iE/P; E = \E\e ±irt> (22) 

and the two fluids have oppositely directed, non-zero components of momentum 
perpendicular to the propagation direction, whose magnitude is related to the electric 
field strength: 

\P±\ 2 = ( 23 ) 
= 7 2 - P || 2 -l (24) 

where \p±\ and \E\ are constants. This extra degree of freedom (compared to the 
subluminal striped wind) brings with it a restriction on the plasma density in the form 
of a non-linear dispersion relation: 

u = 7em^ p (25) 

where u v = (8irne 2 /m) , and 7em = P/ a//3 2 — 1 is the Lorentz factor that corresponds 
to a velocity of c/ /3. Seen from a frame moving radially outwards with this speed, the 
wave is spatially homogeneous and the magnetic field vanishes. 



4-1. Perturbation expansion 

Applying a short wavelength approximation to the above equations reveals that the 
slow radial evolution of the wave amplitude is governed by the conservation of particle 
number, energy and radial momentum. 

Cj 
C F 

Cq 



f 2 (J) 
f 2 (F) 
mc 2 

r 2 (Q) 
mc 



f 2 2cnp\\ 
■ f 2 2nc7P|| + 



2ncp\\ + 



\E\ 2 \ 
4:7imc(3 J 
\E\ 2 / 1_ N 
liimc \ B 2 . 



(26) 
(27) 

(28) 
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The dispersion relation (}2"o"j) implies that the proper density is proportional to 7^. 
Therefore, because of the conservation of particle flux, pu oc 7 EM /? s2 - The constant of 
proportionality is not important for the analysis of the wave properties, but it can be 
found from the spin-down power: 

P|| = ^% (29) 

where is the strength parameter defined in (151). As a result of (1291) . this wave 

propagates only at a sufficiently large distance from the star. Since p\\ < 7 < /x, 

the wave propagates only for 

^l7em ol 

r > > — 

H fi 

= 4.3 x lOVr 1 ^ 2 (4vr/a) 1/2 (30) 

where z/4 = /i/10 4 . In the following, the dimensionless radius will be expressed in terms 
of this minimum radius for propagation, using the variable 

R = rfi/a L (31) 

From ( 1261) and ( 1271) the fx parameter is 

(F) 



mc 2 (J) 



\E\ 2 



= 7 + 5 ' ' 2 (32) 

= 7 + %^ (33) 
PP\\ 

where the dispersion relation ( |25|) has been used in the last step. Similarly, from (|26l) 
and PS]) one finds 

^ = „ + ^( 1 + l) (34) 



The radial dependence of the wave follows from ( 1331) and (IMjl . together with 



Pll = ^ (35) 
7 2 = l+PH 2 +p ± 2 (36) 

/? = -7=^= (37) 
V7em - 1 

Using (134]) to eliminate in (133]) enables one to find a quadratic equation for py, so 
that for a given wave speed 7em, there are at most two distinct physical solutions for p\\ 
and, hence, for the radius. In order to relate these solutions to the striped wind solutions 
that carry the same particle, energy and momentum fluxes, it suffices to express Cq / Cj 
in terms of the magnetisation parameter a using ( fl3l) The radial dependence of these 
two modes for /1 = 10 4 and a = 10 2 is shown in Fig. [TJ 
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Figure 1. The properties of circularly polarised EM waves for fj, = 10 4 , a = 10 2 . The 
quantity p w is the radial momentum of the electron and positron fluids seen from the 
wave frame. It is positive (wave slower than particles) at R = 1, but becomes negative 
at large R. 



To find the behaviour at large radius, first write an expression for p±, using (|34|) 
and (E51): 



P± 



EM 



R 2 Mm " 1) 



Mi 



7em 
R 2 



l + a 
2fi 



(3f 



As — > oo, 7 EM cannot rise faster than R, since pn < \i. Therefore, for both modes, 
p±_ ~ -R -1 , i.e., the wave amplitude decreases as 1/R. 



4-2. Free- escape mode 

Writing 7 EM = yR and expanding in powers of 1/R gives 



1 - 



a 



2fi 



1/2 



as R — y oo 



Assuming /x/(l + a) ^> 1, this mode has the properties 



7EM -> -R/ Vl + cr 



,V2 



(39) 

(40) 
(41) 
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PH -> — — (42) 

as i? — )■ oo. At large distances, its Lorentz factor increases linearly with R, but the radial 
momentum of the particles tends to a constant, so that the density falls of as R~ 2 . The 
wave essentially detaches itself from the particles, and becomes a vacuum wave with a 
frequency large compared to the local plasma frequency. Most of the energy flux at large 
radius is carried by the wave, in the sense that the ratio of the two terms on the RHS 
of (133)) is 7emP± 2 / (Pp\\) ll — °~ ^ 1- Thus, the ratio of electromagnetic to particle 
energy flux tends to the same value as in the MHD wave that carries the same total 
particle, energy and momentum fluxes. 

At smaller radius, this mode changes character, roughly at the point where the 
wave Lorentz factor equals that of the particles, i.e., R ~ ji/ \/l + a. Inside this point, 
the wave moves outwards more slowly than the particles. Since pj « 7 « //, most of 
the energy flux in this inner region is carried by the particles. The wave amplitude p± 
is nearly constant. 



4-3. Confined mode 

The second mode is very different. At large radius, 7em tends to a constant, p\\ — > 
and 7 — > 1. In this case, combining ( 133)) and ( 134")) for large R leads to 




A*-^ = ^-l)J U + (43) 

This equation has one superluminal ((3 > 1) root provided \i > (1 + cr)/2, as is expected 
under pulsar conditions. For large \i it corresponds to 

7em -> W8) 1/4 (44) 
P± 11/R (45) 
PH -> ^ 3/2 / (R 2 V8) (46) 

Close to R = 1, this wave moves more slowly than the particles, but quickly 
overtakes them. However, it does not detach itself. The radial momentum of the 
particles decreases as R~ 2 , so that the proper density remains constant. The ratio of 
the wave frequency to the local plasma frequency approaches the value (/i/8) 1 ^ 4 . Most 
of the energy flux is carried by the wave — for this mode at large radius the ratio of 
the terms in (1331) equals fi, substantially greater than the ratio in the equivalent MHD 
wave. 



5. Summary 

An analysis is presented of some simple nonlinear waves that propagate in the inner 
and outer parts of a pulsar wind. These are described by the equations of relativistic 
MHD in the inner part, and by a model with two relativistic cold fluids (electron and 
positron) in the outer part. The inner solution tends to convert energy flux carried by 
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the fields into particle energy flux as it accelerates outwards. Two electromagnetic outer 
solutions were found, both of which do the opposite, converting particle energy flux to 
wave energy flux. One of the outer solutions tends to a vacuum wave at large radius 
and has the same ratio of particle energy flux to flux carried by the fields as does the 
corresponding MHD wave. The other remains a plasma wave at large radius, but carries 
an even larger fraction of its energy flux in the fields. 

The most important property of the outer wave modes is likely to be their damping 
rates. Previous analyses of the linearly polarised plane wave case have found strong 
damping by radiation reaction for pulsar parameters [23], and, for general parameters, 
strong growth of the Weibel instability driven by the counter-streaming fluids [24) . These 
and similar effects are likely to determine the structure of the "termination shock" at 
the interface between the pulsar wind and the nebula, but have yet to be analysed in a 
global model. 
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